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Corona virus disease (COVID-19) have become a world health problem 
because they have attacked many people worldwide. Because this virus has 
spread massively in almost all countries, including Indonesia, the Indonesian 
government made some policies and rules to close down the hotels and 
restaurants to avoid the spread of COVID-19. Because of that, estimation of 
the number of closed down restaurants and hotels in Jakarta is vital for 
avoiding COVID-19 spreads further to other people, either domestic or 
foreign. In this paper, the adaptive neuro-fuzzy inference system (ANFIS) is 
chosen as the estimation method. In estimating the number of closed 
restaurants and hotels using ANFIS, supporting variables such as the amount 
of casualties in Jakarta, the amount of casualties in Indonesia, and the 
amount of casualties in the world is required. As a result, ANFIS can 
estimate the amount of closed down restaurants and hotels approaching the 
target. The simulations are organized by partitioning the dataset into two 
parts: data of (80%) and data of testing (20%). According to ANFIS 
simulations, ANFIS can estimate the number of closed down restaurants and 
hotels in training data with optimal RMSE equals 0.5324 and testing data 
with optimal RMSE equals 5.3198. 
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1. INTRODUCTION 


Tourism is a vital sector for increasing the revenue of a country. People often do travel around the 
world for vacation, business, or other purposes. Before 2020, the number of incoming tourist worldwide were 
about 532,953 million in 1995, and then in 2019, they can reach 1,461 billion. International tourist arrivals 
have enormous economic impacts. In 2019, contribution of the sector on international tourism is 
approximately USD 8.9 trillion to world gross domestic product (GDP) or 0.3 per cent of the global GDP [1]. 

In the early march of 2020, the corona virus disease (COVID-19) has started to spread in Indonesia 
and become a world health problem because they have attacked many people globally, and some of them 
have died. This disease is firstly found in China and then spread worldwide, including in Indonesia. 
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COVID-19 is the disease caused by the coronavirus resembling severe acute respiratory syndrome (SARS) so 
that it is named SARS-CoV2. The symptoms are fever, tiredness, dry cough, sore throat, difficulty in 
breathing or shortness of breath. The transmission occurs in droplets through sneeze or cough from an 
infected person [2]. All countries implement lockdown, ban foreign visitors, and travel restrictions from other 
countries to reduce the spread of COVID-19. The lockdown causes some people to get reduced earning and 
lost their jobs. Banning foreign visitors and travel restrictions makes hotels’ and lodging revenue decrease, 
and employers can dismiss their workers and employees [3]. 

The impact of COVID-19 causes the growth of economy worldwide to decrease 3% in 2020. The 
growth of economy in developed countries declined by 6.1%. According to the forecast, the world trade will 
decrease between 12.5% and 31.5% in 2020. The global inflation rate is expected to be 3% in 2020 [1]. 

Because this virus has spread massively in many countries, the Indonesian government has made 
some rules to close down the restaurants and hotels to prevent the spread of COVID-19. Jakarta is the city in 
Indonesia where the first victim comes and the city with the most COVID-19 positive victims. Based on 
Statistics Indonesia (BPS) data, the number of hotels and lodging in Indonesia are 28,243. The COVID-19 
impacts the hotel and restaurant so that many of them are closed down about 49.54% [4]. Because of that, 
estimation of the amount of closed down restaurants and hotels in Jakarta is essential for preventing 
COVID-19 from spreads further to other people, either domestic or foreign. In this research, the adaptive 
neuro-fuzzy inference system (ANFIS) will be the estimation method. 

ANHIS leverages a supervised learning algorithm and Takagi-Sugeno fuzzy inference systems. One 
characteristic of an adaptive network is that it consist of some adaptive nodes. In general, each node has 
particular purposes, and the output is influenced by the parameters and signals. In ANFIS, a hybrid method is 
used to train premise and consequent parameters [5]. Both ANFIS and Backpropagation are almost similar 
methods in training and testing the data. They work forward and backwards in training the data for producing 
minimal mean square error (MSE) in making predictions and using maximal accuracy in making 
classification [6], [7]. Furthermore, some optimization methods such as particle swarm optimization (PSO) 
[8], artificial bee colony (ABC) [9], or cuckoo search [10] are used to find parameters in both ANFIS and 
backpropagation. 

In the preceding research, several forecasting methods for estimation have been applied statistically 
by using exponential smoothing method [11], autoregressive integrated moving average (ARIMA) for stock 
price estimation [12], H-infinity for estimating world crude oil price [13], Kalman filter on determining 
trajectory, steam temperature and water level estimation [14], stock price estimation [15], blood supply 
estimation [16]—[18] and fuzzy logic by using ANFIS [19]. Furthermore, backpropagation has been used in 
weather prediction affected by air temperature and humidity [20] and precipitation [21]. This algorithm is 
used in testing and training data with a particular proportion. In this paper, ANFIS estimates the number of 
closed down restaurants and hotels in Jakarta. In calculating the amount of closed restaurants and hotels 
using ANFIS, it requires supporting variables such as the amount of positively infected COVID-19 in 
Jakarta, the amount of dead people caused by COVID-19 in Jakarta, the amount of positively infected 
COVID-19 in Indonesia, the amount of dead people caused by COVID-19 in Indonesia, the amount of 
positively infected COVID-19 in the world, and the amount of dead people caused by COVID-19 in the 
world. According to simulation results, ANFIS can estimate the amount of closed down restaurants and 
hotels in data of training with optimal root mean square error (RMSE) equals 0.5324 and data of testing with 
optimal RMSE equals 5.3198. We also repeat these simulations with different membership functions. 


2. RESEARCH METHOD 

The approach used in this paper is by ANFIS. ANFIS uses a supervised learning algorithm and 
Takagi-Sugeno fuzzy inference system. Furthermore, ANFIS applies a hybrid method in training the data. 
Jang [5] proposed a hybrid method to train consequent and premise parameters. There are two sections of the 
hybrid method, i.e. forward and backward path [22]. 


2.1. Algorithm of adaptive neuro-fuzzy inference system (ANFIS) 
Given two inputs x y and a single output f. Suppose the following rules are used [5]: 

Rule 1: If c is P, and dis Q,, then g, =a,c+fPidt+y, 
Rule 2: If cis P, and dis Qz then gz = a,c + Bod + yz 
with P,, Pz and Q,, Q2 are the membership functions. Parameters @, 61, y; and @2, Bz, yz denote consequent 
parameters. 

The diagram of ANFIS model is shown in Figure 1. Each layer of ANFIS has a different calculation 
given in (1) until (5). In the first layer, each node is the function parameter. The output of every node denotes 
the degree of membership value from the membership functions. The membership function is a Gaussian 
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function u(y) = exp (- (=) ) or a generalized bell membership p(y) = . The parameters a, 8, y 


1 
a) 
of the membership function are called the premise parameters. In the second layer, each node is the strength 
of firing from the corresponding rule. The output node denotes the outcome of signal multiplications, that is 
sent to the following node as in (2). In the third layer, each node is the proportion between the i-th strength of 
firing and the total strength of firing. This outcome is defined as the normalized strength of firing as in (3). In 
the fourth layer, every node is the function defined as in (4). Finally, in the fifth layer, this node aggregates 
whole signals of the preceding nodes as in (5). 
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Figure 1. The diagram of ANFIS model 
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2.2. Hybrid method 

Jang [5] proposed a hybrid method for training of consequent and premise parameters. There are 
two sections of the hybrid method, i.e. forward and backward path. In the forward direction, signals move up 
to layer 4, and then the consequent parameters are obtained using the least square estimation (LSE). The rate 
of error when moving backwards in the backward path, and the premise parameters are obtained using 
gradient descent. A single level of the hybrid method is defined as an epoch. A hybrid method combining 
LSE and gradient descent can result in a faster convergence rate than the Backpropagation method [5]. 

Premise parameters a, B,y are trained and optimized using backpropagation and gradient descent: 
Suppose the network has L layers, and the k -th layer consists of #(k) nodes, such that node Jin the k -th 
layer is (k,i) where the output is OF. Suppose the amount of training data is p, such that the sum of square 
error E,, is in (6): 


a ig 0hs) (6) 
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with T,,, is the target of the m -th component in the p -th data and One denotes the output of the m -th 
component of the layer L by the p -th data. All computations for N training data are obtained using: 


N 
E= > 2 
p=1 


After E, is obtained, calculate gradient descent in the L -th layer in (7). Generally, the node (k, i) is 
calculated using the chain rule as in (8): 


OE 
gf, = ~2(Tup — Op (7) 


k4+1 

OEy ee pele ae dOmp (8) 
ke = +1° kK 
aok, m=1 aon dok, 


foreach 1 <i< #(k) wih1<k<L-1. 

If ~is a premise parameter of the network, then the chain rule is defined in (9) with S is the set of 
nodes in which the outputs influenced by a. All error value for N training data are calculated using (10). And 
update premise parameters using (11) and (12): 


2p = 5 <5 20% ® 
oa = Dea (10) 
Aa = -= (11) 
Gnext = Seurr — AQ = Curr — (— =) (12) 
with 7 denotes the rate of learning. The RMSE is calculated using: 
RMSE = Jere => Sut 6 tee Ang (13) 


with T,; denotes the target and Y,, denotes the outputs. 


3. RESULTS AND DISCUSSION 

Dataset is taken from the report of closed down hotels and restaurants during March 1, 2020, until 
April 30, 2020 (61 days), issued by Indonesian restaurants and hotels association, with the observed area in 
Jakarta as the first city where COVID-19 are first discovered. These data will be estimated by the ANFIS. In 
addition, some data about the number of COVID-19 victims (positive and dead) in Jakarta from March 1, 
2020, until April 30, 2020, is obtained from the official COVID-19 Jakarta website. The number of 
COVID-19 victims (positive and dead) in Indonesia from March 1, 2020, until April 30, 2020, is obtained 
from the Ministry of Health, the Republic of Indonesia. The data of the number of COVID-19 victims 
(positive and dead) worldwide from March 1, 2020, until April 30, 2020, is obtained from the worldometer 
website. They will be used as inputs in ANFIS, with the output is the number of closed down restaurants and 
hotels in Jakarta from March 1, 2020, until April 30, 2020 (61 days) as estimations. 


3.1. Correlation test 

The graph of the number of closed down restaurants and hotels from March 1, 2020, until April 30, 
2020 (61 days) in Jakarta can be seen in Figure 2, where the first closed down hotels occurred on March 17, 
2020, then they are followed by other hotels because the spread of COVID-19 is enormous. After that, the 
highest number of closed down hotels and restaurants occurred from April 1, 2020, until April 5, 2020. The 
next time, the graph shows a minor downtrend because some hotels and restaurants are opened again with 
strict health protocols. 
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Figure 2. The number of closed hotels and restaurants 


For estimating the number of closed down hotels and restaurants using ANFIS as output (y), we 
need supporting variables such as the amount of positive casualties in Jakarta (x,), the amount of dead 
casualties in Jakarta (x2), the amount of positive casualties in Indonesia (x3), the amount of dead casualties in 
Indonesia (x4), the amount of positive casualties in the world (xs), and the amount of dead casualties in the 
world (x6). The graph of COVID-19 victims in Jakarta from March 1, 2020, until April 30, 2020, can be seen 
in Figure 3. There are two parts from the graph with different colours, such as the number of positive victims 
and dead victims. From the chart, the amount of COVID-19 positive victims and the amount of dead victims 
in Jakarta from March 1, 2020, until April 30, 2020, increase exponentially. For example, the number of 
positive victims in Jakarta can reach 4138 on April 30, 2020, and the number of dead victims in Jakarta can 
get 381 on April 30, 2020. 


The number of Covid19 victims in Jakarta 


the number of victim 


Figure 3. The number of COVID-19 victims in Jakarta 


The graph of COVID-19 casualties in Indonesia from March 1, 2020, until April 30, 2020, can be 
seen in Figure 4. There are two parts from the graph with different colours, such as the number of positive 
victims and dead victims. From the graph, the amount of COVID-19 positive victims and the amount of dead 
victims in Indonesia from March 1, 2020, until April 30, 2020, increase exponentially. The number of 
positive victims in Indonesia can reach 10,118 on April 30, 2020, and the number of dead victims in 
Indonesia can get 792 on April 30, 2020. 
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Figure 4. The number of COVID-19 victims in Indonesia 


The graph of the amount of COVID-19 victims in the world from March 1, 2020, until April 30, 
2020, can be seen in Figure 5. There are two parts from the graph with different colours, such as the number 
of positive victims and dead victims. From the diagram, the amount of COVID-19 positive victims and the 
amount of dead victims worldwide from March 1, 2020, until April 30, 2020, increase exponentially. For 
example, the number of positive victims in the world can reach 3,299,603 on April 30, 2020, and the amount 
of dead victims in the world can get 233824 on April 30, 2020. 
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Figure 5. The number of COVID-19 victims in the world 


For determining linear relationship among supporting variables such as the amount of positive 
casualties in Jakarta (x7), the amount of dead casualties in Jakarta (x2), the amount of positive casualties in 
Indonesia (x3), the amount of dead casualties in Indonesia (x4), the amount of positive casualties in the world 
(xs), the amount of dead casualties in the world (6) against the amount of closed down restaurants and hotels 
in Jakarta as output (y), we compute the correlation in Table 1. After that, we use the linearity test using the 
t-test to determine the linear relationship [23] where the t-test is in (14): 


= rvn-2 (14) 


t 


with n is the amount of data and r denotes the correlation coefficient and the t-distribution having n — 2 
degrees of freedom. 
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Table 1. The correlation between some inputs against the number of closed down restaurants and hotels 


Correlation variables Correlation value t-test p-value Linear relationship 
Positive (Jakarta) The closed hotels and restaurants 0.7525 8.7764 0.0000 Yes 
Dead (Jakarta) The closed hotels and restaurants 0.7631 9.0685 0.0000 Yes 
Positive (Indonesia) The closed hotels and restaurants 0.7127 7.8050 0.0000 Yes 
Dead (Indonesia) The closed hotels and restaurants 0.7303 8.2114 0.0000 Yes 
Positive (world) The closed hotels and restaurants 0.8084 10.4165 0.0000 Yes 
Dead (world) The closed hotels and restaurants 0.7633 9.0765 0.0000 Yes 


From Table 1, all the amount of positive casualties in Jakarta (x), the amount of dead casualties in Jakarta 
(x2), the amount of positive casualties in Indonesia (x3), the amount of dead casualties in Indonesia (x4), the 
amount of positive casualties in the world (xs), the amount of dead casualties in the world (xs) have linearity 
relationship against the amount of closed down hotels and restaurants in Jakarta as output (y). 


3.2. Estimation results 

Before applying the estimation method, we require to divide the data into two parts: training data 
and testing data. The proportion is 80% training data and 20% testing data [24]. Figure 6 displays the division 
between testing data (red plus marks) and training data. We make this pattern for training and testing data 
because the data is nonstationary. Data used in training and testing data must have resembling characteristics 
and properties [25]. 


the number of closed hotel 


time 


Figure 6. Data divided into two parts training data (circle) and testing data (plus) 


The parameters in the simulation of ANFIS are: 


— Maximum of epoch : 30 
— The amount of membership functions : 2 in each input 
— The amount of rules : 64 


— The amount of linear parameters (consequents) : 448 
— The amount of nonlinear parameters (premises) : 36 

The simulation results of ANFIS are shown in Figures 7 to 10. Two membership functions are used: 
the circle curve denotes ‘small’ and the green curve denotes ‘large’. Figure 7 displays the optimized 
membership functions of inputs such as the amount of positive casualties in Jakarta, the amount of dead 
casualties in Jakarta, the amount of positive casualties in Indonesia, the amount of dead casualties in 
Indonesia, the amount of positive casualties in the world, and the amount of dead casualties in the world. 
Each membership function leverages a generalized bell (gbell) containing three nonlinear parameters. 
Figure 8 displays the convergence process of RMSE. In the early epoch, and the RMSE was pretty big. 
During the optimization process, the RMSE is decreased to 0.5324. Figure 9 displays estimation results (the 
number of closed down restaurants and hotels) on training data. It seems that there are the comparison and 
error between the target and the output. From the training process until the 30 epoch, it is obtained optimal 
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ANHIS parameters. Then, optimal ANFIS parameters are used in testing data. Figure 10 displays estimation 
results (the number of closed down restaurants and hotels) on testing data. 
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Figure 10. Estimation result on testing data of the amount of closed down hotels and restaurants in Jakarta 


From the simulation in Figure 9 and Figure 10, we achieve the estimation with the RMSE in (13) are: 
— Data for training : 0.5324 

— Data for testing : 5.3198 

They seem that error on the testing data is more significant than on the training data. 
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We replay the simulations five times using other membership functions such as gaussmf resembling 
normal or gaussian distribution shape, trapmf using trapezoidal shape function, trimf using triangular shape 
function, pimf, and dsigmf. Finally, the results are displayed in Table 2, ie. RMSE on ANFIS with different 
membership function with maximum epoch is 30. The best membership function for training data is dsigmf 
with an RMSE equals 0.4422, while the best membership function for testing data is Trapmf with an RMSE 
equals 3.0560. They seem that errors on the testing data are more significant than on the training data in all 
five membership functions. 


Table 2. The RMSE on ANFIS with other membership functions 


Membership function Max epoch RMSE on training data RMSE on testing data 
Gaussmf 30 0.5375 4.8326 
Trapmf 30 2.1666 3.0560 
Trimf 30 1.8291 3.9111 
Pimf 30 2.1729 8.4132 
Dsigmf 30 0.4422 6.4076 


4. CONCLUSION 

In estimating the number of closed down restaurants and hotels using ANFIS, it is required 
supporting variables such as the amount of positive casualties in Jakarta, the amount of dead casualties in 
Jakarta, the amount of positive casualties in Indonesia, the amount of dead casualties in Indonesia, the 
amount of positive casualties in the world, and the amount of dead casualties in the world. ANFIS can 
estimate the number of closed down restaurants and hotels approaching the target. Simulations are conducted 
by dividing the dataset into two parts: training data (80%) and testing data (20%). From ANFIS 
computatitons, ANFIS can estimate the amount of closed down restaurants and hotels in data of training with 
optimal RMSE equals 0.5324 and data of testing with optimal RMSE equals 5.3198. The improvements of 
this research are making estimation and classification the results by data mining and machine learning 
techniques. 
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